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Field ofl the Invention 

Tkc present imvcnLian relarttos to a process for increasing the efficiency of a 
computer in finite; alemiont simulations by efficient automatic construction ofl 
suitable basis functions fan computation of approximate solutions, and to a 
5 computer for performing tWat process. 

Background ofl t hie Invention 

A plurality of technical and physical! phenomena can: be described by partial 
differential equations. They include among others problems from fluid mechanics 
1Q (for example, flow ground an airfoil}, electromagnetic fi'eld theory (for example. 
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electrical field behavior in a transistor)' or elasticity theory (for example, 
deformation of a car body). Accurate ktoowilodgc and description of sucli processes 
are a central element in the construction and optimization of technical objects. Tb 
save time-consuming and cost-intensive experiments, there is great interest in 
5 computer-aided! simulations. Finite element processes (FE processes) haw become 
established and have been tjhle topic of intense research for a long time. This also 
applies to automatic mesh generation processes as a foundation for construction of 
suitable basis functions. 

1Q Figure 1 illustrates in the left half the prior art in the process of FE simulation for 
a linear boundary value problem as a typical model example. Proceeding from the 
data describing the geometry of the technical object to be simulated, first a system 
of baisis functions is constructed which on the one hand enables fulfilment of 
boundary) conditions and on the other is suited for approximation of the unknown 

15 solution. Then, using these basis functions a system of! linear equations is set up 
by numerical integration methods. Finally, the coefficients of the unknown 
approximation are determined as the solution of this system of linear equations. 

FIE methods ar thoir use axe the subject matter of a series of patents. For example, 
20 U.S. Patent Ka. 4, 81i9 5 161 discloses a system where FE approximations of a large 
class of differential equations are automated. U.S. Patent N!o. 5,731,817 discloses 
a process fox generation of hoxahedral meshes forming the foundation for a FE 
simulation process. 

25 In most FE processes of practical relevance basis functions are used which are 
defined on a decomposition! produced by generating a mesh of the simulation 
region. Figure 2a shows a selection of conventional elements; their dimension, 
degree, smoothness and panametors are listed in Figure 2b. A survey of meshing 
methods of planar regions can he found for example in K. Ho-Le, Finite Element 

30 Mesh Generation Methods: A nervicwi and classification. Com. Aided Design 20 
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(111988), 27 - 38. Generating aj mesh' far complicated three-dimensional regions is 
extremelyi difficult using the current state oil knowledge, as shown by S. Qwcn, A 
survey of unistruotured Mesh Generation Technology;, Proceedings, 7th 
International Meshing Round Table, Sandia National Dab (1998J, 239 - 2S7L The 
5 processes require extensive amounts of computer time and are to some extent not 
yet oompletely automated. But recently, there has been a series of very; innovative 
new approaches. For example, AJ. Euchs, Optimierte Delaunayi-Triangulierungen 
zur Vcrnetzung gptrimmter NIURJBS-Korper, University of Stuttgart, 1999, 
simulates a force distribution in order to achieve an optimum distribution of 

IQ triangulation points. In U.S. Patent No. 5,729,670 two- and three-dimensional 
meshos are produced by, solving flow problems; this is an interesting reversal 1 of the 
conventional FE mechanism. In addition, many algorithms have been developed to 
improve individual aspects of mesh generation processes. For example DE 
196 21 434 Al and U.S. Patent No. 5, 774, 696 describe a process for elimination of 

15 intersections with prescribed edges or boundary surfaces in Delaunay 
triaiigulations. 

Meatless FE-micthods to date have not acquired any importance for applications. 
Rothi in the Lagrange multiplier method (see for example J.H. Bramble, The 
20 Lagrange Multiplier Method fbr Dirichlet's Problem, Math. Comp. 

37(|1981), 1 - 11), and also in the penalty method (see for example P. Bochev and 
M. Gunzburgcr, Finite Element Methods ofi Least Squares Type, SI AM review 
40 (11998], 789 - 837), tihe treatment of boundary conditions represents a major 
problem in the use af simple, stable basis functions. 

25 

Tn many technical simulations, automatic mesh! generation is very complex and 
requires by far the largest p;art of the computer time. Furthermore, the 
approximation power of the conventionally used linear and multilinear basis 
functions is low. To achieve accurate results, a large number of basis functions 
30 must he used, and thus, a cmrraspondingly large system of equations must be 
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solved. Higher order trial functions on triangulations generally likewise have an 
unfavorable ratio between the attainable accuracy and the number of basis 
functions used. Finally, smooth basia functions cannot be easily defined on 
unstructured meshes. Very special constructions are necessary] already for 
5 continuously differentiate elements (see Figure 2a). 

Summary of thie Invention 

The abject of the present invention is to increase the efficiency of known FE 
methods md computers which carry out FE methods by efficient construction of 
basis functions with favorable properties. In particular the meshing of the 

1Q simulation region will be completely eliminated, optional boundary conditions are 
fulfilled, accurate solutions are obtained with relatively few coefficients, and the 
resulting system of equations will he solvable efficiently. In this way the 
disadvantages of the prior art will be overcome, and thus, the accuracy and speed 
of the simulation of physical properties in the engineering and optimization of 

15 technical objects will be improved. 

Some contrail terms and notations which are used in the following description of 
the process as claimed in the invention will be explained first. 

20 The simulation region Q is a bounded set of dimension d = 2 or d = 3 on which the 
physical quantities to be studied will be approximated by means of FE methods. 
The boundary of the simulation region is denoted by F. A grid with grid width h 
is defined as a decomposition of a subset of the plane or the space in grid cells Z k . 
Depending on the dimension d each grikl cell is a square or a cube with edge length 

25 h. More precisely, Z% = kh 4- [IQ. k] d , where M belongs to a set of integer d-vectors. 
The uniform tensor product B-spllilnes in. d variables of degree n with 1 grid width h 
are denoted fty sec for example 0. de Door, A' Practical Guide to Splines, 
Springer, 1978. They arc functions which can be continuously differentiated (n — 11) 
times and which an the grid cells agree with polynomials of degree n 5 as shown in 
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Figures 4a and 4b. Pilgpirc 4ai sham, flic support of f he Rl-splinc of degree 
n = 2, dimension d = 2 andl smoothness m — 1. In Figure 4h, the resulting tensor 
product] B-spline ^ is shown. The support i.e. the union, of all grid cells, on 
which the B-spline bk is not identical to zero, consists of (n + 1 ] d grid cells; more 
5 precisely, Q h = kH-ti [0, (n + l)h] d . In all figures, the B-spline 6 fc is marked ait the 
point kh y i.e. for example in the case d = 2 at the lower left corner of the support. 
Fbn EE simulations only those B-splines are important which have support 
intersecting the simulation region 0; they are called relevant B-splines. The 
relevant B-splines are again divided into two groups^ those B-splines for which the 
lid part of the support inside the simulation region is larger than a prescribed bound si 
are called inner B-splines. All other relevant B-splines are called outer B-splines. 

Tlhc object of the present invention is achieved by the process defined in claim 1. 
Special embodiments of the invention are defined in the dependent claims. Claim 
15 11 defines a computer system as claimed in the invention. 

The right half of Figure 1 shows the incorporation of the process of the present 
invention into the counse of a FE simulation iin the prior art and the substitution 
of certain process steps of a FB simulation, in the prior art by the process of the 
20 present invention. 

Tnput 1 of the simulation region Q can be done via input devices, in particular also 
by storage of data derived from comiputer-aidod engineering (0AD/|CAM). For 
example, the data used! in the engineering of a motor vehicle can be incorporated 
25 directly into the FE simulation 3n the present invention. 

In- input 2 and storage mflthe type of boundary conditions natural and essential 
boundary conditions arc distinguished. The basis in the present invention is 
constructed for homogeneous boundary conditions of the same type. In particular, 
30! for essential boundary conditions, the basis functions vanish on the boundary R 



Inihomogeneous boundary cowditions can he treated in tlic assembly of tlie F1F 
systems using methods which correspond to the prior art. 

Finally the control parameters are read in 3. They; relate to the degree n and the 
5 grid width h of the B-splines to be used and the bound s for classification of the 
inner and outer H-splines. If specifications are omitted, all these input parameters 
can be automatically determined by evaluation of merit functions which arc 
constructed empirically or analytically. 

1Q The following construction of the basis functions of the invention is divided into 
the steps shown schematically in Figure 3, which will now be described. 

After neading in the simulation region S7, in the first process step a grid covering 
the simulation region Q is generated. Then it is checked which of the grid cells lie 

15 entirely inside, partially] inside or not inside the simulation region fl The cell 
types 4 are determined, and tihis information about the cell types is stored. This 
essentially requites inside/outside tests and determinations of intersections 
hctwiecn the boundary F of the simulation region Q and the segments or squares 
which bound the grid ccllls. Figure 7 shows the input and output data for this 

20 process step. 

Tn the second process step, using the information about the cell types, the relevant 
B-splines arc first determined. Then 1 the classification 5 into outer B-splines is 
performed; the corresponding lists of indices are denoted by / and J\. To this end, 
25 tlic size of those parts of the supports ofl the B-splines which lie within the 

simulation region is determined using the data obtained in the first process step, 
and compared wiithi the prescribed Hound a. Figure 9 shows the input and output 
data far this process step. 

30 In. the third process step, coupling coefficients etj are computed 6; they join the 



inner and outafl B-spllines according to the rale 

Bi(x) = bM + J} *iA * e i. dij 

j€J(i) 

Hence, an extended B-spline Bj is assigned to each immer R-spline Bj. The 
construction and the properties of the index sets ,T{i\ and the coupling coefficients 
5 e» j are given as follows. The index sets J\{i) consist of indices of outer E-splines. 
They; correspond to complementary index sets of indices of inner. B-spiimcsj 
i.e., i belongs to Ifyj) if and only if j belongs to J(i). For a given outer index j the 
index set is an array , i.e., a quadratic or cubical arrangement of (n + l) d inner 
indices which is characterized by, a minimum distance to the index j. For, a given 

10 outer index j and an inner index i in the index set J(j), let Pi be the d-variate 
polynomial ofl degree n in each variable which has the value 1 at the point i and at 
al other points of the array the value 0. Then thle coupling coefficient j is 
given as the value of pi at the point j\ y i.e., ay = Pi{j\ The specific values of the 
coupling coefficient can- either be tabulated for different degrees and relative 

15 positions of j and or can be easily computed using Lagrange polynomials. 
Figure 11 shows the input and output data for this process step. 

Tf natural boundary conditions are given, the extended splines defined in equation 
(1) arc used without 1 modification for further implementation of the FB process. M, 
20 on the other hand, essential boundary conditions are given, a weighting according 
to the rule 

BiW^^mix), iei (2) 

has to be performed. The pertinent interrogation, takes place in an optional process 
step 6a. The functions defined in this way are called weighted extended B-splines 
25 OWEB-splincs). Formally, the extended B-splines used under natural boundary 
conditions correspond to the special! case w(x) = 1. Th'ey are, therefore, also called 
WEB-spImcs. For the case of essential boundary conditions, the weight function w 
is characterized as follows! Pan all points <r. of the simulation region, w(x] can be 
bounded from above and below, by posiltilve constants, which are independent ofl a;, 
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times the distance dist(#) of the point x from the boundary T\ Tn other words, w 
is positive within Q and tends to zero in the vicinity of the boundary T as fast as 
the distance function dist. For simulation! regions wihJiah arc bounded by 
elementary, geometrical objects (circles, planes, allipses. etc.) a suitable weight 
5 ftinction can optionally be given in explicit analytic for.mi. Otherwise, eomjpuitaiticm 
rules should be used which typically represent a smoothing of! the distance 
function. The scaling fttctor Ifwfei) is calculated by evaluating the weight 
function at the weight point aa$. This can be any point in the support of the 
D-spline flfe which is at] least half! the bound a/2 away from the boundary. 

Id 

As a result of the prooess of the present invention, a computation rule for the 
WEB-spIines Bi (compare definitions (1) and (2)) is obtained which have all 
favorable properties according to the objectives of the invention. Thus the EE 
method can be continued according to the prior art. But, in doing so, it is possible 
15 to advantageously exploit the regular grid structure of the basis functions 
according to claim 11. 

In summary, for the prooess of the present invention the coupling of outer and 
inner B-spIine is, among other things, important. As a result s the constructed 

2Q basis has the properties which aire essential for FE computations. In particular, a 
basis Bi (ji from the index set I) according to the present invention, is stable, 
uniformly with' respect to the grilcl width /?,, and the error has the same order as for 
tlic R-spiincs b& in the approximation of smooth functions which satisfy the same 
boundary conditions. On the other hand, the fulfillment of essential boundary 

25 conditions is ensured by using the weight function w. 

Other objects, advantages amd salient features of the present invention will become 
apparent from the followring detailed description which, taken in conjunction with 
the drawings, listed below, discloses preferred embodiments of; tih'e present 
3d invention. The features mentioned in the claims and in the specification can be 
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essential to th!c invcntiioii individually or in any combination. 



Brief Description of th!e Drawings 

Referring to the drawings which farm a part of this disclosure] 

5 Figure 1 is .a flow chart of the individual steps of the prior, art in the course 
of the firuite element simulation and integrates the determination of 
the W1EB basis into this process; 

Figure 2a compares certain finite elements of tihe prior art to the WEB 
element shown in Figure 2b and lists tihe parameters relevant to 
HQ finite element approximations; 

Figure 3 is a flow chart of the process steps for determining the WIEB basis; 

Figures 4a audi 4b show a support and the corresponding tensor product 
B-spIine of degree 2; 

Figures 5a and 5b illustrates the problem formulation of the first 
15 embodiment (displacement of a membrane under constant 

pressure) and of the corresponding solution; 

Figure 6 shows the cell types for tihe first embodiment; 

Figure 7. surveys the input and output data ofl the process for determining 
the cell types] 

20 Figure 8 illustrates, by way ofl example, the classification of the B-splines 
for the first! embodiment); 

Figure 9 outlines the input and output data of the prooess for classifying 
the B-splines; 

Figure 10 shows the coupling coefficients of an outer B-spline and the 
25 corresponding tamer B-splincs for the first embodiment; 

Figure 11 surveys the input and output data of the process for computing 
the coupling coefficients; 
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Filgures 12a and 12b illustrates the construction of tWc weight function of 
the preferred embodiment; 

Figure 13 shows the support of a; WER-spline and the corresponding 
coupling coefficients for the first emibodimient; 

5 Figures 14a and 14b explains the problem formulation of a second 

embodiment (incompressible flow) and its solution using the flaw 
lines and the distribution of the flow velocity]; 

Figures 15a to 15c show the changing classification of B-splines for the 

B-spline degrees n = 1, 2, 3 and the same grid width for the second 
liQ embodiment; 

Figures 18a to 116c provide information about the error development in the 
finite element approximation using WEB-splines and about the 
computing time behavior of WEB approximations for the second 
embodiment) 

15 Figures 17a amd 117b compare the WEB basis to a prooess based on linear 
trial funiotJons on a triangulation (prior ant); and 

Figure 1& shows a computer system according to the present invention. 



Detailed Description of the Invention 

One especially favorable embodiment of the process of th'c present invention, called 
the WEB process, is determined by the following specifications. 

9 The bound s is chosen such that the inner B-splincs are characterized by, requiring 
that at least one of the grid cells of their support lies completely in the simulation 
region Q. Since to determine the relevant B-splines, the intersection of the grid 
cells and the boundary] T must be computed anyway, the classification requires no 
significant additional computing time. The weight point] x$ is chosen as the 
lid midpoint of a grid cell in the support of the B-spline hi which lies completely in 
the simulation region SI. This is also efficiently possible since the determination of 
one such: cell is already part of the classification routine. 

Tf no explicit, analytic form ofl the weight function is known, it is defined by 



Figures 12a and 12b illustrate the construction of the weight function. Here, the 
parameter S indicates the width of the strip Sl§ within which the weight ftonction 
varies between the value 0 ot the boundary ofl the simulation region and the value 
1 on the plateau omi Sl\Q^ u The parameter & is chosen such that the smoothness of 
2Q the weight functrmn is ensured. 

One important advantage of the prooess is that no meshing ofl tihe simulation 
region is necessary. In technical aipplicatilons, this results in clear savings of 
computing time and storage capacity and simplifies the course of thle simulation. 
25 The process structure for two- and three-dimensional problems is formally and 
technically largely identical. This enables time- and cost-saving implementations of 
solutions for diverse applications based on uniform program structures. The use of 
Bi-splincs corresponds to the industrial standard in the modeling of geometrical 
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objects, and thus forms a natural connection between FE and GAD/|GAM 
applications. Extensive existing program libraries from both! fields can be used for 
implementing a EE simulation based un thJe process of the present invention. The 
basis functions constructed using the WEB process have all standard properties of 
5 finite elements. This includes especially the stability of the basis. It implies, for 
example, that for linear elliptic boundary value problems the condition! number of 
the resulting system of equations does not grow faster than for optimal 
triangulations as the grid width becomes smaller. For applicationSj this means, for 
example, that linear systems of! equations as they typically arise in FE methods, 

10t can be efficiently solved by iterative algorithms. Furthermore, for a given degree, 
the appnoxilmaition order is maximal and the number of necessary parameters 
minimal. Thus, very accurate approximations are possible with a relatively small 
number of parameters. Specifically, this can mean that the accuracies which so far 
required the use of mainframe computers cam now be achieved with workstations. 

15 The regular grid structure of the basis of the present invention permits a very 
efficient implementation, especially for assembling and solving FE systems. 
Moreover, by using the weight function, the boundary conditions can be satisfied 
during simulation without affecting the regular grid structure of the basis 
functions. Finally* by using multigrid methods to solve the linear systems arising 

20 in linear elliptic boundary value problems, one can achieve that the overall solution 
time is proportion al to the number of unknown coefficients, and thus optimal. 

The process of the present invention in the special preferred embodiment (WiEB 
process} is illustrated using the first embodiment shown in Figures 5a and 9b. The 
25 differential eqiuation and boundary conditions are chosen to be very elementary so 
that in addition to the construction of the WEB basis of the present invention the 
entire course of the FE simulation can be fallowed without major additional effort. 

Figure 5a shows an clastic membrane, which is fixed along the edge F of a planar 
3Q simulation region Q., and on which a constant pressunc / = T acts inside the region. 
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Wlitih suitable norraaJftaaliam, the displacement u satisfies the Pbisson equation 
wiiithl homogeneous boundaryi conditions, 

—Aw = 1 in Q 

r ii = 0 onE 

o 

The displacement u or the deflection of the membrane is depicted im Figure 5b. As 
described above, the WIEB process is divided into the following steps. 

Input 1 of the simulation region fib The boundary Y is a periodic spline curve of 
liQ degree 6, which is stored by| its control points 20 (in Figure 3a identified with black 
dots). 

Tnp.ut 2 of the boundary conditions: the homogeneous boundary condition is 
essential so that the construction of a weight function is necessary. 

15 

Tnput 3 of tho control parameters: The degree n — % and, in order to make the 
figures easier, to understand, a relatively large grid width h = 11/3 are used. 

Determination 4 of the cell types: As illustrated iln Figure 6, the simulation region 
20 is covered! by a gnid 21, which contains the grid cells of the supports of all 

BI-spKncs of potential relevance fbr the basis construction. The type determination 
in the example yields 69 outer grid colls 22, and 11 inner grid cells 24, and 20 grid 
colls 23 on the boundary. 

Classification S of the B-splinesi Here, the support of the B-spline 6*. is a square 
Qk with corners 

AfcJA, (jfci + 3 5 h)h, {h +. 3, & -H 3)/i, (k lf k 2 + 3) A"; 

25 ^(-^o) and 0(2,1) arc shown- in Figure 8. The grid points kh of the relevant 

rS-splincs, far which intersects the Interior of the simulation region, are marked 
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in Figure 8 by a point or a circle. Aft grid points ih for inner R-splincs (ji from thic 
index list J), for which at least one cell of the support Qi lies entirely within Q, are 
manoked by, a point. For example, far i = (-4, 0). the grid cell (-2, Ojft +" [0, /i] 2 lies 
entirely in flL All grid points jh for outer H-splines (j from the index list] J), far 
5 which no cell of the support Q$ lies entirely in are marked by a circle. 

Computation 6 of the coupling coefficients: To determine the coupling coefficients 
eij fior each fixed j of thle index list J the nearest 3 x 3-array 

= {fti, 4 + Mi + 2)] x {4, £2 -H 1, 4 + 2} 

Id of indices in I is sought. In Figune 10, for the outer grid point j = (-11, 2), which is 
marked with a circle, the array is identified with points. Figure 10 likewise 
shows tihe corresponding coupling coefficients in a matrix representation. They are 
computed Ky hivardate interpolation. For example, the interpolating polynomial: 
for i = (-1, -1) is 

15 Pi(x) = -x 1 (x 1 -H 2)x 2 (x 2 - l)/2. 

Its value at the paint x = j = (—1, 2) is Pi(ji) = I = e^. One notices that many, of 
the coupling coefficients are 0. Thlis is a typical phenomenon. The coupling 
coefficients aire not equal to 0 for alii of the index list only if the indices i 
arc different from tlhc index j in each component. 

20 

Computation rule 7 for the weight functiiom The weight function is given by 
equation (|3)| with n = 2 and S = 0.2. The parameter & is computed numerically. It 
must be small enough to prevent singularities of the distance function. To compute 
the distance function, the process generates a program which uses Newton's 
25 method. Since the weight function is not equal to 1 only in a boundary strip, the 
complexity in the sublscqjuanit evaluations is low. 

Output: Figure 13 shows the support of a WRB-spline B { and the data necessary 
far its description. These are the indox list Jfli) ofl thle outer B-splines bj coupled 
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with! Bi> the coupling coefficients Cy, and the weight' point These data arc used 
in conjunction with the weight funnt'ion for generating' the computation rule for 
the WEB-splines. 



5 The further course of the FE simulation follows the prior art. 



Assembly 9 of the FE system: The entries of the system matrix amd of. the 
right-hand side are 



a. 



si = / gnadE fc gradBi, F k = / fB ky 
Jn Jn 



k. i e K 



liQ Tlhe system af equations QO = F for the basis coefficients d in this example has 
dimension 31. The matrix entries are computed using numerical integration , 
likewise the integrals F k > 



Solution 10 of the FE system: The Galerkin system is solved iteratively with the 
15 conjugate gradient method with SSOR preconditioning used to accelerate 
convergence. After 24 iteration steps the solution is found within machine 
accuracy (jtoJai-anoe <| le — 14). 

Computation 11 and output 12 of the approximation: The approximation 
20 computed with the process as claimed in the invention ils u = Jjj an d is 
shawm graphically im Figure 5b. The relative error of the L 2 -novm is 0.028. 

The efficiency of the process of the present invention in the special preferred 
embodiment (jWEB process) is illustriated in a second embodiment using the 
25 simulation of an incompressible flow. The arrangement of two circular obstacles 
shown in Fligurc 14a in a channel! wilth parallel boundaries serves to illustrate the 
principal strategy. For complicated geometries, as are typically present in specific 
applications, the process works completely analogously and efficiently. In 
Figure 14a the stream lines 25 ame shown within the region bounded by I\ to T 4 - 
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and by V& and F 6 . The differential equation is: 

Au - 0 3n Q. 

with the boundary conditions 

du _ du ^ du ^ ^ n 

— = u 0 on Td , — = - V0 Qr\r 2 , — = 0 om r 3s . , . , P 6 * 
an 9n an 

o Thle flow velocity v — — gradu is shown in tilie bottom halfi of thle figure. 

The construction of the WBB basis of the present invention proceeds completely 
analogously to the first embodiment. The sole difference is that a weight function 
is mot necessary because of the natural boundary conditions. 

10. 

Figures 15a to 115c show the classification of the relevant B-splines for different 
degrees n (sec also Figure 8). In the figure, the inner B-splines ftj, which are taken 
into the WEB basis without extension, are marked by solid triangles. For small h 
the number of these B-splines increases, i.e., Bi = hi for most of! the WEB basis. 
15 For degree = 3 this is the case for 236 of 252 indices in the example. 



Figure 16a shows in two diagrams the numerically determined relative I^-error. ofl 
the potential (left half of the figure) as a function of the grid width hi = 2~ k with! 
k = 1, ... , 5 and the numiemcalliyi estimated order ofl convergence m (right halfl ofl 
20 the figure]. Hero, for different degrees of the WEB-spline the fallowing markers are 

used: * (!rc - *)> 0 (™ = 2 ): A (fa = 3), □ (n = 4) and * (n = 5). As expected, 
m fts n + 1, i.e., m approximate error reduction by a factor 2 ra+1 when the grid 
width is cut in half. Analogously, for the relative approximation error of the flow 
velocity shown in Figure 16b (jff^-norcm of the solution, left Half ofl the figure), an 
25 order of convergence m^n flrighfi half of figure) is obtained with an associated 
error reduction tiy roughly ai factor 2 n when the grid width is out in halfL 

Figure 16c (right half of the figure) shows the computing time in seconds for 
construction' of the WEB tiasis as a function of the number of resulting basis 
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functions, measured on a Plenitiuim TT processor wiithl 400 MHIz. For example, for 
construction of a WEB basis of degree 3 with' grid width A = 0.125 with 2726 
WEB-splines 1.32 seconds are necessary. One notices that the complexity for 
generating the WEB basis is largely independent of the degree n of the basis. In 
5 the left half of Figure 16c the number of OG-itcraitions relative to thto number of 
basis functions is shown. Thus, for the corresponding system with 2726 unknowns, 
65 POG-iterations are required. The total computing tilme including assembling 
and solving the Gialerkin system is roughly 2.48 seoonds. 

10 Figures 17a and 17b compare the WEB process with a standard solution process 
which meshes or triangulates the simulation region (Fig. 117a) and uses hat 
functions. The graph shows in Fig. 17b the L2-error relative to the number of 
parameters. The results of. the standard solver are marked with boldfaced 
diamonds and are compared to the results achieved using the WEB basis of 

15 degrees U to 5. For examp1e ? an accuracy of 10~ 2 is achieved with the WEB 

process by using 21 3 basis functions with degree 2 and an overall computer time of 
0.6 seconds. To achieve the same accuracy^ the standard method with linear hat 
functions required Q657 basis functions. 

20 In the assessment of the standard solution process two other aspects must be 
considered. On the other hand, Figure 17b illustrates that even a moderate 
accuracy of 10~ 3 cami only be achieved with hat functions when far more than one 
million coefficients are used. This shlowis that when using hat functions accurate 
results generally require an enormous computing and storage capacity or cannot be 

25 achieved at all with the prior art. On the other hand, the complexity required for 
meshing increases with' the complexity of the simulation region. Tn contrast to 
realistic applications, the region studied Here is comparatively; simple to 
triangulate due to its simple structure. 

30 The two-dimensional example shows the performance gain by the WEB process. 
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An even greater increase in performance is passible in three-dimensional problems. 
On the one hand, the complexity for mieshing, whicli is eliminated in the WEB 
process, is much greater. On the other haind. the reduction in the number of 
required basis functions becomes much' more noticeable than in the 
5 two-dimensional case. 

Eilgure 18 shows a device according to the present invention, especially a computer 
system 30, with input devices 311, 32 3 33, output devices 34 and a control unit 35 
which controls the course of! the process. To carry out the process of the present 

1Q invention and in particular for purposes of parallelizatian of. the pertinent 

computations, the central control unit 35 preferably uses several arithmetic logic 
units dALU) or even several central processing units (CPU) 36. These allow 
especially parallel processing for, the process steps classification 5 of the B-splines, 
in particular also intersection of the regular grid with the simulation region Q, 

15 determination 6 of the coupling coefficients ey s and/or evaluation of the weight 
function w(x] at points x of the simulation region £1 

The computer units 36 here access the common data resources of the storage unit 
371. The data aain be input, for example, by a keyboard 31, a machine-readable 

20 data medium 38 via a corresponding read station 32 and/|or via a wire or wireless 
data network with a receiver station 33. Via the read station 32 or a pertinent 
data medium 38, the control program, which controls the process execution, can be 
input, and, for example, can bo permanently filed on the storage media 37. 
Accordingly, the output devices 34 can be a printer, a monitor, a write station for 

25 a machinc-rcadahlc data medium amd/on the transmitting station of a wire or 
wireless data network. 

While various embodiments have been chosen to illustrate the invention, it will be 
understood by those skilled in the art that various changes and modifications can 
30 be made therein without departing from the scope of the invention as defined in 



